Defect Fluctuations and Lifetimes in Disordered Yukawa Systems 
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We examine the time dependent defect fluctuations and lifetimes for a bidisperse disordered 
assembly of Yukawa particles. At high temperatures, the noise spectrum of fluctuations is white and 
the coordination number lifetimes have a stretched exponential distribution. At lower temperatures, 
the system dynamically freezes, the defect fluctuations exhibit a 1// spectrum, and there is a power 
law distribution of the coordination number lifetimes. Our results indicate that topological defect 
fluctuations may be a useful way to characterize disordered systems. 
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B 



X3 
O 

> 
O 

in 
o 
\o 
o 



I 

o 
o 



X 



I. INTRODUCTION 

In two dimensional systems of interacting repul- 
sive particles such as dusty plasmasii^ and colloidal 
assemblies^'^'^, heterogeneous particle motions have 
been observed in the dense liquid phase. The hetero- 
geneity appears close to the disordering transition where 
the system changes from having predominantly hexago- 
nal ordering to being heavily defected. Recently, it was 
demonstrated that in the disordered regime where dy- 
namical heterogeneities are present, the time dependent 
density fluctuations of the topological defects exhibits a 
1// power spectrum^. The heterogeneous motion dis- 
appears at higher temperatures, and simultaneously the 
noise spectrum changes to a white form and the noise 
power is strongly reduced, indicating a lack of correla- 
tions in the fluctuations. It was argued that motion in 
the regime with dynamical heterogeneities occurs in a 
correlated manner and is spatially concentrated in areas 
containing a higher density of topological defects. As a 
result, the creation and annihilation of topological defects 
is strongly correlated, leading to the 1// noise signature. 
Recent experiments on two-dimensional colloidal systems 
exhibiting heterogeneous motion have confirmed that the 
colloids with more rapidly changing coordination number 
are correlated with the regions of motion^. 

In the systems mentioned so far, the ground state at 
low temperatures is an ordered hexagonal lattice; how- 
ever, it is known that intrinsically disordered or glassy 
systems also exhibit dynamical heterogeneities^iS,. A nat- 
ural question to ask is how do topological defect fluctu- 
ations behave in systems that are inherently disordered 
and that have no transition to an ordered state as the 
temperature is lowered. Although there is no ordering 
transition in an intrinsically disordered system, there can 
be a temperature at which a dynamical slowing down oc- 
curs, giving rise to glassy type behaviors. It would be in- 
teresting to study how or if the defect fluctuations change 
as dynamical freezing is approached. In the disordered 
system, the length scale of the dynamical heterogeneities 
increases as the temperature is lowered'^'^. If the moving 
regions are associated with more highly defected regions, 
this may also be reflected in the defect fluctuations. 



In this work we examine a two-dimensional system of a 
bidisperse mixture of Yukawa particles. Previously, this 
system has been shown to form a disordered state at all 
temperaturesiS. We analyze the time dependent defect 
fluctuations, noise spectra, noise power, and the coordi- 
nation number lifetimes for varied temperatures. At high 
temperatures where there is no heterogeneous motion, 
the defect fluctuations have a white power spectrum, the 
coordination number lifetimes are very short, and there is 
an exponential decay in the distribution of lifetimes. As 
the temperature is lowered, the defect fluctuations show 
a 1// power spectrum with increased noise power. We 
also find a peak in the noise power at a finite tempera- 
ture we label T„. Near this temperature, the motion is 
highly heterogeneous and the coordination lifetimes are 
power law distributed. For temperatures below T„, the 
noise power is reduced and the coordination number life- 
times become extremely long. We also probe the system 
by examining the motion of a single particle driven by an 
external force. At T = there is a well defined threshold 
force for motion^fi. At finite temperatures, creep occurs; 
however, there is a well defined kink in the velocity force 
curves at a threshold force which vanishes at T„. These 
results suggest that measuring topological defect fiuctua- 
tions may be a useful probe for understanding glassy and 
jamming behaviors in disordered systems. 



II. SIMULATION 

We consider a two-dimensional system with sides of 
length L and with periodic boundary conditions in the 
X and y directions. The system contains N particles in- 
teracting via a Yukawa or screened Coulomb potential. 
For any two particles i and j of charge qi and qj located 
at positions and r^ , the pair interaction potential is 
= qiQj ex.p{-Kr,j)/rij. Here rtj = jr^ - r^l and 
1/k is the screening length which is set equal to 2.0. 
We consider a 50:50 mixture of particles with charges 
qi and 52 where (71/(72 = 1/2. We have previously 
shown that this system forms a disordered assembly at 
all temperaturesiS. The particles obey overdamped dy- 
namics and the time dependent particle positions and 
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FIG. 1: (a) Voronoi construction for a snapshot of a bidis- 
perse particle system at T = 1.5. White polygons have 6 
neighbors, while black and gray polygons have 5 and 7 neigh- 
bors respectively, (b) The particle positions (black dots) and 
trajectories (black lines) for the same system during 10^ sim- 
ulation time steps. 



velocities are obtained by integrating the overdamped 
Brownian dynamicsii equations of motion, which for a 
single particle i is 



dt 



(1) 



Here rj is the damping coefficient and = 
— Vil^(ry). The particle-particle interactions are 
cut off at lengths larger than 2/k and are calculated 
using a neighbor index method for computational effi- 
ciency. Further cutoff lengths have negligible effects. The 
thermal force is modeled as random Langevin kicks with 
(Ff ) = and (F'^(i)F^(t')) = 2i^kBT5{t - t'). The ini- 
tial conditions are obtained by simulated annealing. For 
our parameters, if all the charges are set equal to the 
larger charge ^2, melting as defined by the onset of de- 
fects occurs at a temperature of T = 4.0. The driving 
force Fd is applied to a single particle only, such that 
Fd = /dx for the driven particle and F^ = for all other 
particles. For most of this work, we set fd — 0. We 
anneal from a temperature of T = 8.0 down to a final 
temperature in steps of AT — 0.25 and sit at each incre- 
ment temperature for over 10^ Brownian dynamics (BD) 
time steps. After annealing, we analyze the topological 
defect density and particle motions at a fixed tempera- 
ture. Unless otherwise noted, we consider systems of size 
L — 32 containing N — 830 particles. 

In Fig. 1(a) we show the topological defects as obtained 
from a Voronoi construction for a system at T = 1.5. 
The Voronoi construction forms a polygon around each 
particle and is used to define the coordination number 
Cn of the particle. Particles with five, six, and seven 
neighbors corresponding to Cn = 5, 6, and 7 are marked 
black, white, and gray, respectively. The system is also 
highly defected at other temperatures and there is little 
change in the average defect density with temperature; 
however, the time dependent fluctuation rate of the de- 
fects changes significantly with temperature as we show 
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FIG. 2: (a) Power spectra S{i',T) obtained from the time 
series of the topological defect density fluctuations Pe{t) ob- 
tained at r = 1.5 (upper curve) with a fit (solid line) of 
a — 0.98 and T = 4.5 (lower curve) with a fit (solid line) of 
a = 0.1. The T = 4.5 curve has been shifted down for presen- 
tation purposes. Inset: the fitted values of a obtained from 
the power spectra vs T. (b) The noise power So vs T. Inset; 
So for T — 1.0 to r = 2.5 for varied system sizes L — 16 
(squares), 32 (diamonds), 48 (triangles), and 54 (circles). 



later. In Fig. 1(b) we illustrate the particle trajectories 
during 10^ BD time steps at T = 1.5. The motion is 
highly heterogeneous, and certain regions show motion 
during this time interval while other regions do not. At 
higher temperatures the motion becomes more homoge- 
neous throughout the sample. 



III. DEFECT FLUCTUATIONS AND POWER 
SPECTRA 

In order to characterize the defect fluctuations, we per- 
form a series of simulations at different T and examine 
the time dependent density of six-fold coordinated par- 
ticles, Pe{t), defined as Peit) = YlfLi SiC'Jt) - 6), 
where is the coordination number of particle i. The 
power spectrum of the resulting time series is defined as 



Peit)e 



-2T:iut 



dt 



(2) 



The noise power Sq is defined as the average value of the 
noise spectrum at a particular value of the frequency u. 
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After annealing we wait 10^ time steps before acquiring 
the time series Pe(t) to avoid any transient behaviors. In 
Fig. 2(a) we show the power spectra of Peit) for T = 1.5 
and T = 4.5. At T = 1.5 the power spectrum has a 
1//" form with a ~ 0.98, as indicated by the upper soHd 
hue. For T — 4.5 the noise spectrum is approximately 
white and a « 0.1. We plot the evolution of a with T 
in the inset of Fig. 2(a). For T < 2.0, we find 1// noise 
with a ~ I, while at high T, the noise spectrum becomes 
white with a w 0. The particle trajectories also indicate 
that the motion is no longer heterogeneous for T > 2.0. 
In monodisperse systems near the disorder transition, a 
similar 1// noise signal was observed when the particle 
motions were highly heterogeneous^. We note that for 
T < 1.5 there is very little particle motion and the defect 
fluctuations are strongly suppressed. This can be more 
clearly seen by examining the noise power Sq at fixed 
= 5 as shown in Fig. 2(b). At high temperatures, 
So is low, but it increases rapidly to a peak value as 
the temperature is reduced to T = Tn = 1.5, and then 
decreases quickly for T < 1.5. We associate the peak 
in Sq at T„ with the temperature at which the system 
begins to undergo a dynamical freezing. 



In the inset of Fig. 2(b), we plot the variation in 5'o 
for systems of different sizes ranging from L = 16 to 
54. As L increases, the peak in the noise power becomes 
more pronounced. Additionally, although the height of 
the peak in Sq increases with L, at higher or lower T 
the noise power at the non-peak temperatures is reduced 
with increasing L. For temperatures away from T = 1.5, 
the individual power spectrum curves have a rollover to 
a white spectrum at the lowest frequencies. At T = 1.5, 
the a = 1.0 behavior persists down to the lowest fre- 
quencies, as shown in Fig. 2(a). The suppression of the 
noise power away from T = 1.5 in the larger systems can 
be understood as arising from the fact that there is a 
length scale I associated with the correlated creation or 
annihilation of defects. The fluctuations with the largest 
I correspond to the lowest frequencies. If I is constant, 
then when the system size is small or comparable to I, 
the fluctuations will have a 1// character down to the 
lowest observable frequencies and the noise power will be 
increased away from the peak in Sq, as seen for L = 16 in 
the inset of Fig. 2(b). As the system is made larger, the 
fluctuations begin to average out, a low frequency cutoff 
appears in the 1// spectrum, and the noise power is re- 
duced. Near the dynamical freezing temperature T„, the 
fluctuations grow with the system size and 1 // noise is 
present for all arbitrarily large L and low j/. For T < T„, 
the fluctuations are again cut off and the noise power is 
reduced. We have compared our results against simula- 
tions of collections of non- interacting particles, where the 
motion is homogeneous and diffusive at all temperatures. 
In this case, there is no peak in 5*0 and the spectra are 
white for all temperatures. 
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FIG. 3: (a) The measure of the particle displacements d{t) vs 
t for (top curve to bottom curve) T = 3, 2.5, 1.5, 1.0, and 0.5. 
The dashed line is a linear fit showing the behavior expected 
for a system undergoing normal diffusion, (b) The value of 
d{t) at t — 10^ time steps measured for the same system as a 
function of T. 



IV. DIFFUSION 

To quantify the amount of motion present in the sys- 
tem as a function of temperature, we measure the diffu- 
sion of the particles by computing the square displace- 
ments of the particles from an initial position as a func- 
tion of time, 

JV 

d{t)=J2\h{t)-r,{to)\\ (3) 

i=l 

where ri{t) is the position of particle i at time t with 
the effect of crossing the periodic boundary conditions 
removed. For particles undergoing normal diffusion, 
d{t) oc t. In Fig. 3(a) we plot d{t) at T = 3.0, 2.5, 1.5, 
1.0, and 0.5 from top to bottom. The data was taken 
over a period of 10^ time steps which started after the 
system had entered a stationary state. In general, we find 
a transient time period which lasts for 10^ time steps af- 
ter annealing. This transient behavior was excluded from 
the analysis so that > 10^ time steps. For T > 1.5, 
the long time behavior of the system is consistent with 
normal diffusion, as indicated by the dashed line showing 
a linear slope in Fig. 3(a). For T < 1.5, within our time 
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FIG. 4: (a) The coordination number C„ vs time for a single 
particle at T = 1.5 showing occasional changes, (b) Cn vs 
time for a single particle at T = 3.5 showing much more 
rapid changes, (c) The histogram of the coordination number 
lifetimes P{ti) for T = 3.5 (lower curve) and T = 1.5 (upper 
curve). The solid line is a power law fit with /3 = 1.6. Inset: 
P{ti) at T = 5.0, 3.5, 2.5, and 1.5, from left to right. 




FIG. 5: The velocity 14 vs fd for a single driven probe particle 
for T = 3.5, 2.0, 1.25, 1.0, and from left to right. The peak in 
the noise power of the defect fluctuation power spectra occurs 
at T = r„ = 1.5, which also corresponds to the onset of the 
subdiffusive motion. Inset: the threshold force /c extracted 
from the velocity-force curves vs T. 



frame d{t) saturates to a constant value, while at T = 1.5 
the long time behavior is more consistent with subdiffu- 
sion, d{t) oc t" where a < 1.0. We note that this agrees 
with the measures of the noise fluctuations which indi- 
cate that the 1// noise which appears at T = T„ = 1.5 is 
associated with the presence of system- wide slow dynam- 
ics. For T < 1.5 the fluctuations are instead associated 
with the rattling motions of single particles rather than 
long time system wide motions. In Fig. 3(b) we plot the 
value of d obtained at the end of 10^ steps. For T < 1.5, d 
saturates to a constant value since within our time frame 
the particles are no longer diffusing significantly, while d 
increases with T for higher T. 

We note that since we are restricted to performing sim- 
ulations over finite periods of time, we cannot be certain 
that the subdiffusive motion at T = 1.5 does not cross 
over to regular diffusion at extremely long times. If this 
were the case, the power spectrum of the corresponding 
very long time scries would show a low frequency cutoff in 
the 1// signature at the time scale where the regular dif- 
fusion begins to dominate. What we can conclude from 
our data is that over the time regimes which we arc able 
to access, the appearance of the 1// defect fluctuations 
is correlated with the subdiffusive behavior. 



V. DEFECT LIFETIMES 

We next consider the lifetimes of the coordination 
numbers. In Fig. 4(a) we show the coordination number 

Cji as a function of time for one particle in a system at 
T = 1.5. During this time frame, the particle is predom- 



inantly sixfold coordinated with C„ = 6, while there arc 
occasional jumps to a configuration with C„ = 5 or 7. In 
some cases, such as for t > 7000, the particle maintains 
Cn ^ Q for a longer period of time. We obtain a similar 
plot regardless of which particle in the system is selected. 
In Fig. 4(b) we plot C„ vs time for the same particle at 
T = 3.5, showing that the particle coordination number 
is changing much more rapidly. The coordination life- 
times ti are obtained by measuring the length of time 
that elapses between changes in C„ for a given particle. 
We histogram this quantity over all the particles, and plot 
the resulting P{ti) in Fig. 4(c) for two systems at T = 3.5 
(lower curve) and T = 1.5 (upper curve). The straight 
solid line is a power law fit to P{ti) oc with /3 = 1.6, 
while the curve at T = 3.5 is fit well to a stretched expo- 
nential. In the inset of Fig. 4(c) we plot P{ti) for T = 5.0, 
3.5, 2.5, and 1.5, which shows that for the higher temper- 
atures the average lifetime {ti) is shorter. For T = 5.0, 
P{ti) can be fit reasonably well to an exponential distri- 
bution. For the intermediate temperatures, we can fit the 
data with a stretched exponential form, and at T = 1.5 
we obtain a power law fit. For T < 1.5, the lifetimes be- 
come extremely long and we observe close to a bimodal 
distribution where there are a few particles jumping back 
and forth giving rise to short values of ti, while C„ of the 
other particles does not change for the duration of the 
simulations over times of 2.5 x lO"" BD steps. 
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VI. YIELDING BEHAVIOR 

The results from Figs. 2, 3 and 4 suggest that near 
T — 1.5, some type of dynamical crossover or freezing 
occurs. We can also probe the system by slowly driving 
a single particle through it. If the system acts rigidly at 
low drives, then there is a threshold force fc above which 
the driven particle moves with respect to the background. 
Previous work on this system at T = has character- 
ized this threshold forc&lS.. Experiments have also shown 
evidence for a threshold force to motion for disordered 
colloidal systems^. In general, we find that creep can 
occur at low but finite temperatures; however, a thresh- 
old force can still be identified by a sudden increase or 
kink in the particle velocity vs applied force. In Fig. 5 we 
plot the velocity vs applied drive fd for a single probe 
particle which has a charge q/ qi = 5.0 moving through 
a disordered assembly of particles at temperatures rang- 
ing from T = 0toT = 3.5. If the surrounding par- 
ticles are absent, the probe particle moves linearly with 
V (X Fd- The velocity- force curves are obtained by slowly 
increasing the applied drive fd- For T < 1.5, 14 shows 
a sharp downward concavity at finite fd- We define the 
threshold force fc as the value at which = 5 x 10~^. 
For T > 1.5, the curves show a slight positive concav- 
ity and if we linearly extrapolate the curves to = 0, 
the threshold force is indistinguishable from zero. In the 
inset of Fig. 5 we plot the threshold force fc extracted 
from the velocity-force curves. We have considered var- 
ious rates of increasing fd and find no differences in the 
extracted values of fc- The threshold decreases with in- 
creasing temperature up to T = 1.5, and above this tem- 
perature there is no threshold for motion. These results 
add further evidence that there is a dynamical freezing 
or jamming that occurs near T = 1.5, as also reflected 
in the defect fluctuations. We note that it is beyond the 
scope of this paper to determine the true nature of the 
finite temperature dynamical freezing, such as whether it 
resembles a glass transition. 

VII. DISCUSSION 

The main goal of this paper is to introduce a new mea- 
sure, the fluctuations in the defect density, for character- 
izing disordered systems near freezing. The question of 
whether there is really a finite temperature freezing tran- 
sition in the equilibrium system near T = 1.5 is beyond 
the scope of our paper. Although we have found several 
signatures for some form of dynamical slowing down at 
T — 1.5, it is possible that for very long times beyond 
our simulation time scales, the system would show linear 
diffusion and the fluctuations would lose the 1/f char- 
acteristic at very low frequencies. We note that there 
is other evidence that topological fluctuations and de- 



fects appear to be correlated with jamming. Recently, 
the shape features of topological defects have been uti- 
lized to study the jamming transition in two-dimensional 
granular materials^^. In this case the fluctuations are 
nonthermal; however, a jamming transition can still oc- 
cur. Measurements of fluctuations of Pe{t) for the gran- 
ular system, in the manner we suggest in the present 
work, reveal the same trends that we have found here. 
The defect fluctuations have a white noise characteristic 
away from the jamming transition which crosses over to a 
1/f signature near jamming with a peak in the finite fre- 
quency noise power—. This indicates that the topologi- 
cal fluctuations near jamming or freezing may universally 
exhibit 1/f noise features. 

There has also been recent work on two-dimensional bi- 
nary glass forming systems showing correlations between 
the onset of glassy properties and the disappearance of 
certain types of topological defects at finite temperature. 
This adds further evidence that there may indeed be a 
finite temperature freezing or jamming transition in two- 
dimensional disordered systems that is connected with 
the motion of defectsi^. 



VIII. SUMMARY 

To summarize, we have shown that the topological de- 
fect fiuctuations is a useful quantity that can be employed 
to understand dynamical freezing and heterogeneities in 
disordered systems. We specifically find that in regimes 
where there are dynamical heterogeneities, the time de- 
pendent defect density fiuctuations show a 1// charac- 
teristic noise spectrum, and that at higher temperatures 
where the motion is uniform, a white noise spectrum ap- 
pears. The defect noise power also shows a peak just 
before the dynamics of the system freezes, and this ef- 
fect is robust for increasing system sizes. The coordi- 
nation number lifetimes show a power law distribution 
in the heterogeneous regions which crosses over to an 
exponential form at higher temperatures where the het- 
erogeneities are lost. We correlate the peak in the noise 
power with the onset of subdiffusive behavior as well as 
with the temperature where a threshold force for the mo- 
tion of a driven probe particle disappears. Our results 
should be testable for disordered charge-stabilized col- 
loidal assemblies, dusty plasmas, and other disordered 
charged systems. We note that recent experiments on 
granular system near jamming have found the same type 
of behavior in the time series analysis of the defect fiuc- 
tuations as the jamming transition is approached. 

We thank M.B. Hastings for useful discussions. We 
also thank A. Abate for sharing his unpublished data 
with us. This work was supported by the U.S. DoE under 
Contract No. W-7405-ENG-36. 
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